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We study the stability of the nuclear matter in a non-local Nambu-Jona-Lasinio 

' model. We work out the equation of state in a relativistic Faddeev approach where we 

■ take into account the internal structure of nucleon. We show that the binding energy 

C . 

, saturates when the nucleon as a composite particle made of quarks is incorporated. 

' After truncation of the two-body channels to the scalar and axial-vector diquarks, 

(N ■ 

a relativistic Faddeev equation for nucleon bound states in medium is solved in 

(N ; 

^ . the covariant diquark-quark picture. We investigate the nucleon properties in the 



o 
o 

in 
o 



nuclear medium such as the role of diquarks within the nucleon and the in-medium 
modification of the nucleon mass and size. 
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■ I. INTRODUCTION 

Chiral type models such as the NJL model have been successful in explaining the low 
energy physics of mesons and nucleons in vacuum utilizing the concept of the spontaneous 
chiral symmetry breaking Application of the model to low density, however, has posed 
serious problems like the absence of saturation This is in contrast to non- 

chiral models like the Walecka model Q which has been phenomenologically successful in 
describing the nuclear matter ground state in the mean field approximation. 

In chiral models with a Mexican hat potential, the nuclear medium moves the minimum 
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of the vacuum effective potential towards smaller values. This implies a smaller curvature, 
i.e. the scalar sigma mass decreases. Consequently, the attraction between nucleons due to 
sigma meson exchange increases. This effect may destroy the stability of nuclear matter. 
On the contrary, in non-chiral models e.g. the Walecka type model, where the potential is 
simpler, a stable nuclear matter can be found. 

At the nuclear matter saturation point the Fermi momentum kf and the pion mass 
are of comparable scale {kf ~ 2Q2MeV ~ 2771^^). Therefore, chiral dynamics may 
play an important role in a nuclear matter saturation mechanism. Chiral symmetry is an 
approximate symmetry for quarks which is not necessarily translated into a symmetry for 
nucleon degrees of freedom. Furthermore, the delta-nucleon mass splitting A = 293 MeV 
is also comparable to the Fermi momentum kf at nuclear matter saturation. It has been 
shown that the main mechanism behind the delta-nucleon mass splitting stems from the 
fact that the delta and the nucleon have different internal quark structure^ jliSl- Therefore 
one may naturally be led to ask if the internal structure of nucleon is important in order to 
describe the empirical nuclear matter ground state. Quark degrees of freedom have also been 



important in order to descri^ 
GeV and the EMC effect M 



oe deep-inelastic scattering at momentum transfers of several 
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Moreover, at the moment a direct calculation of the many nucleons system based on 
QCD itself is not feasible. Therefore, effective quark theories are the best tools to develop a 
microscopic understanding of a single nucleon together with nuclear and quark matter and 
their phase structure in a unified framework. 

It has been shown that a non-local covariant extension of the NJL model inspired by 



the instanton liquid model 12] can lead to quark confinement for acceptable values of the 
parameters [l^- There are other advantages of the non-local version of the model over the 
local NJL model, for example: the constituent quark mass is momentum-dependent, as also 
found in lattice calculation The regulator makes the theory finite to all orders in the 



correct lor is [1,'tI |. This model Jias 
IS andbaryons 12, 3 



loop expansion and leads to small next-to-leading order corrections 1^ 
been phenomenologically very successful to describe mesons 

in vacuum. More applications can be found in Refs. Q,|20|. Here, we employ this model 



^ The Delta is made of only axial vector diquarks while nucleons can be constructed from both the scalar 
and axial vector diquarks. 
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to investigate the role of the internal dynamics of nucleon in a nuclear matter environment. 
We explicitly construct nucleon degrees of freedom by binding diquarks and quarks with 
relativistic Faddeev equations to describe nuclear matter. The relativistic Faddeev approach 

nnnrrnn 

to describe the nucleon as a composite object has been very successful la la llll lia 1211 122, 

HQ- 
Very recently Bentz and Thomas ^ have shown that local NJL models may lead to satu- 
rating nuclear matter when the quark confinement effect is incorporated to avoid unphysical 
thresholds. Although the idea is interesting, there may be some short-coming in their ap- 
proach. Quark confinement in their prescription is incorporated by introducing an infrared 
cutoff in the model. This way quark confinement is permanent and there is no confinement- 
deconfinement transition in their model. In their hybrid approach, the nucleon is described 
as a relativistic bound state in the static approximation^ of a scalar diquark and quark. It 
has been already noticed that the attractive axial-vector diquark plays an important role in 



description of a single nucleon |Sl,llSi- Therefore, in order to treat the internal structure of 
the nucleon and its modification in the medium adequately, we will take into account both 
the scalar and the axial-vector diquark channels and solve the Faddeev equation without 
invoking the static approximation. 

We start by investigating the role of diquarks in the bound nucleon. The in-medium 
modification of the nucleon mass and the nucleon size will be calculated. We obtain the 
nuclear matter equation of state when nucleon is taken as a composite object made of 
diquark and quark. 

This paper is organised as follows: In Sec. II we introduce the model, discuss its pionic 
sector and fix the parameters. In Sec. Ill we set up the description of the nucleon based on 
a diquark-quark picture in the relativistic Faddeev approach. In Sec. IV we consider nuclear 
matter with composite nucleons and present our numerical results. Finally, a summary and 
outlook is given in Sec. VI. 



^ In the static approximation of the Faddeev equation the quark exchange kernel is taken momentum 
independent, which may give rise to an unphysical in- medium singularity due to a reduction of the 
constituent quark mass. 



II. FORMALISM AND SOLUTION IN THE VACUUM AND QUARK MATTER 



We consider a non-local NJL model Lagrangian in terms of quark degrees of freedom with 
SU{2)f X SU{3)c symmetry 

C = ll;[i^ -m^)tlj + ^g^ja{x)ja{x), (l) 

where rric is the current quark mass of the u and d quarks. First, we restrict the interaction 
terms to four-quark interaction vertices in the quark-antiquark channels. We will add in 
section III the diquark channels in order to model the quark sub-structure of the nucleon 
in the medium. The current ja{x) has the scalar {0~^,T = 0) and pseudoscalar (0~,T = 1) 
components 

ja{x) = j d'^xid^xsF^x - X3)F{xi- x)i'{xi)Tai'{x3), (2) 
with Fq, = (l,i75r). The r^'s are flavour SU{2) matrices with tr(rjr,) = 2Sij. Notice that 



there are different approaches j25 



26l | to introduce non-locality into the interactions. We use 



the four-way separability of the non-local interaction which is also present in the instanton 



liquid model 



This assumption considerably simplifies the calculation. Non-locality 



also emerges naturally in the presence of other gluonic field configurations within the QCD 
vacuum <23] and in models based on Schwinger-Dyson resummation techniques |25|. 

In the medium the repulsive vector mesons are important. For simplicity we add a local 
vector field interaction —g^{ip'-)^ipY to the Lagrangian. Locality of this term is permissible 
as long as the interactions in Lagrangian are not fixed by some underlying theory through 
a Fierz transformation. Upon standard bosonization we obtain 



C = J d xid X3ijj{xi) — rric) 6{xi — X3)6{x3 — x) + F{x — X3)F{xi — x) yj{x) + n{x)i^, 

- 'y^u^{x)6{xi - X3)6{x3 - x) ^(xa) - ^[(r^{x) + vr^(x)] + -^uj^{x)uj^{x). 

Lagrangians Q and (jS)) are equivalent in the classical limit. The equations of motion for 
bosonic fields lead to 



(0-(x),7f(x)) = 97r J d'^Xid^X3F{x - X3)F{xi - x)ip{xi) {l,i-f5T)ip{x3), 

u>'{x) = 2g,i^{x)-f^ijix). (4) 

The Lagrangian (Q) can be immediately recovered from the bosonized version Eq. Q by 
substituting Eq. (UJ into Eq. (0). 
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We define the Fourier transform of the form factor by 

F{x-x.) = I ^e-^(— ')-V(p). (5) 

The function f{p) influences the momentum distribution of quarks in the non-perturbative 
vacuum and will be defined later. Next, we perform a mean field approximation by expanding 
the meson fields around their expectation values a = (a), uf = (uj^) neglecting meson 
fluctuations. The mean values of the pion field and space components of the vector field 
vanish in vacuum and for a baryonic matter at rest due to symmetry. The mean field 
Lagrangian in momentum space can then be written as 

2 2 

C = ^ip) {i^ - Mip) - 7%) Hp) - ^ + 7^' (6) 

where the momentum-dependent quark mass and quark propagator are given by 

M{p)=m, + f\p)a, (7) 
5-i(p)=:^-M(p)-7%. (8) 

We treat the quark fields at one-loop level and the meson fields at tree level. This approx- 
imation is consistent with the leading-order behaviour of the model in a l/N^. expansion. 
Within this approximation the effective potential can be calculated: 



where 7g = 2NcNf is the degeneracy factor for quarks. A gap equation is obtained by 

an 

f d^k Mik)pik) 



requiring ^ = 



where the dynamical quark mass is given in Eq. (jZj). 

Our model contains four parameters: the current quark mass m^, the cutoff (A) hidden 
in the form factor f{p) and the coupling constants g^, , g^. The vector field coupling g^ is 
treated as a free parameter and will be adjusted in the medium. We first fix gT^ and rric for 
arbitrary values of A by fitting the pion mass and the pion decay constant to their 
empirical values at zero baryon density. In this way, we can consider the entire parameter 
space of the model. The pion mass is given by the zero of the IPl two-point function Il.^.,^{q) 
for the pion which in random phase approximation has the following form 

n^,(g) = l + 2Tr / ^ f(kh,T,Sikh,T,Siq + k)fiq + k) = 0, (11) 
9n J (^7r)4 
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where q denotes the total momentum of the quark-antiquark pair . The pion decay constant 



is obtained from the couphng of the pion to the axial-vector current 

/. = ^ / 7f^Tr[^75^(^(P-))75r„(5(p+))]/(p-)/(p+) 
J [Ztt)'^ 2 
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9.9.,, f dk ^^^^^^^^ I A^Tr[r'^75/b-)/(p+)5(p-)75r.5(p+)] 



2ml J (27^)4 

X [f\k) + f\p.)) - f{p+)f{p.)f{k) {f{k + q) + f{k - q))l (12) 

where g-^qq is the pion-quark-antiquark coupling constant and related to the corresponding 
loop integral Eq. |TT| by g~qq = {^^^)\g2=mi- We define p± = p ± |g. Notice that due to 
non-locality the one-pion-to-vacuum matrix element gets an additional contribution. This 

n 

extra term is very important to maintain the Gell-Mann-Oakes-Renner relation jl3|] and 
makes a significant numerical contribution. 

In our numerical treatment of the model, we evaluate loop-integrals like that in 
Eqs. ()10I11I12|) in Euclidean space^. We choose for the form factor in Eq. © a Gaussian 
function in Euclidean space, 

f{pE) = exp{-pl/A'), (13) 

where A is the cutoff of the theory. This choice respects Poincare invariance and for certain 
values of the parameters it leads to the quark "confinement" in the sense that the dressed 
quark propagator has no poles at positive p"^ in Minkowski space and consequently, quarks 
do not appear as asymptotic states. Quark " confinement" can occur if the following relation 
is satisfied between the parameters of the model 



a 1 (y/mg + A^ + m,)^ , 

> o exp — . (14) 



V^^^^fTA2 _ 2 " y 2A2 

Notice that the quark propagator has many pairs of complex poles, both for confining 
and non-confining parameter sets which can be conceived as a remnant of the underlying 
confinement. Confinement in this model has direct consequence for the pion, but condition 
Eq. (fT^ is not derived by using any assumption about pions. 



^ The symbol Tr denotes a trace over flavour, colour and Dirac indices. 

^ Notice that here (e.g in Eq. (|10|l ) a standard transcription rules from Minkowski to Euclidean momentum 
(fc" iluji —k^) and reverse has been assumed. So far this assumption has been phenomenologically 

correct ElElElElElEiEiH- 
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FIG. 1: The zero-momentum constituent quark (left panel) and quark condensate (right panel) as 
a function of cutoff for fixed /jr = 93 MeV and 771,^ = 140 MeV when rric 7^ 0. 
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FIG. 2: The constituent quark mass M(p) as a function of momentum p for the three parameter 
sets given in Table 1. 

The emergence of complex poles for quark p rop agat ors was also noticed in Schwinger- 



Dyson equation studies in QED and QCD 



29 



|3l|. It has been recently shown that 



mryons become more compact in the presence of imaginary poles in the quark propagator 
isl ]. In another study it has been noticed that the mass-like singularities are located on 
the time axis if one removes the confining potential in QED in 2+lD, and the mass-like 
sing ularities move from the time axis to complex momenta if there is a confining potential 
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TABLE I: The parameters A, g-,^ and nric for the sets A, B and C, fitted to /t^ = 92.4 MeV and 
m,r = 139.6 MeV. The resulting values of the constituent quark mass M(0) (in brackets we give 
the values for zero current quark mass) and the corresponding quark condensate for the different 
parameter sets are also given. 

Parameter set A set B set C 

A (MeV) 1046.8 1008.1 847.8 

5,r(GeV-2) 3-^_g 352 55.8O 

nic (MeV) 7.9(0) 8.5(0) 11.13(0) 

M(0) (MeV) 297.9(250) 313(264.5) 351.6(300) 

(MgY) 215(208) 212(205) 191(186) 



TABLE IL The first two sets of poles of the quark propagator =bMR it iMj in Minkowski space 
with rric 7^ . 



set A 


set B 


set C 


±391 MeV 


±507 ± i55 MeV 


±408 ± i238 MeV 


±675 MeV 


±1530 ± i253 MeV 


±1575 ± i307 MeV 



We refer to appendix for details of numerical evaluation of loop integrals as of Eq. (fTTjl . 
Solutions of the gap equation ()10|) are shown in Fig. 1. In the left panel the constituent 
mass at zero momentum M(0) is shown as a function of the cutoff. It is obvious that for 
very small cutoff, there is no solution for the gap equation. On the right panel of Fig. 1, 
we show the corresponding values of the quark condensate. These values are within the 
limits extracted from QCD sum rules [s^ and also lattice calculation js^. Notice that the 
trend of the model parameters in Fig. 1 is slightly different to that obtained in the local 
NJL model with a sharp cutoff js^. In contrast to the local NJL model, the dynamical 
quark mass Eq. ((Zj) is momentum dependent (see Fig. 2) and follows a very similar trend to 
lattice simulations Fig. 2 shows that at low virtualities the quark mass is close to the 
constituent mass while at large virtualities it approaches the current mass. 

We analyse three sets of parameters, as indicated in Table IH Set A is called the non- 
confining parameter set since it does not contain imaginary poles of the quark propagator. 
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Sets B and C lead to the quark "confinement" in our convention with complex poles. We 
define the pole mass to be given through the equation 

P^,i, = M\P,,,,), (15) 

where pole solutions are denoted Ppoie = ±M/j±iM/. Equation (fT3j) is a non-linear equation 
and can be solved numerically. The position of the quark poles are given in Table |n] where 
it is seen that for the confining sets the quark poles lie in the complex plane. In contrast to 
the local NJL model, a non-zero current quark mass of order 8.5 MeV leads to an increase 
of about 48.5 MeV in the zero-momentum dynamical quark mass (see Table 1 set B). It is 
also noted that the real part of the pole mass is bigger than the zero-momentum constituent 
quark mass M(0) for both the confining and the non-confining sets, see Tables 1,2. 



III. NUCLEON INTERNAL STRUCTURE 



In order to have a physically realistic description of matter at low density, one has to 
construct nucleons and nucleon matter from the quark degrees of freedom since quarks are 
confined at very low density. In this section, we take into account the quark structure of 
the nucleon in a diquark-quark picture. We first build up a single nucleon out of the quark 
degrees of freedom by solving the relativistic Faddeev equation, then (in the next section) 
we construct nuclear matter in a mean field approximation by means of these individual 
nucleons. In this way, any reference to the quarks will be naturally hidden in the nucleons. 

In the vacuum baryons and diquarks have already been solved in this model by one of 
the authors [18j. Here, we shortly recapitulate the procedure to construct the nucleon as a 
diquark-quark bound state, more details can be found in 18|, |24 1 . 

In order to describe the nucleon as a bound state of a diquark and a quark, we firstly 
introduce diquarks in the model. We truncate the quark-quark interaction to the scalar 
(0"*", T = 0, 3) and the axial vector (1+, T = 1, 3) colour antitriplet. 



Cls,a = gsJs{x)Js{x) + gaJa{x)Ja{x), 



(16) 
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FIG. 3: Top panel: A graphical representation of the Bethe-Salpeter equation for the qq T-matrix in 
RPA approximation. The solid lines denote the dressed quark propagators Eq. (jS)) and double lines 
denote diquark propagators D^^D^'^. The filled circles denote the diquark vertices ^,x^'^- Down 
panel: The coupled effective Bethe-Salpeter equation for the nucleon vertex function (j) is defined 
by amputating the external quark and diquark propagator from the wave function (j) = S^^ D'^tp. 

where the currents are defined 



Js{x) = 


j d'^Xid'^x^f{x 


- X3)fiXi 


-x)ij{x,)[j,CT2P^]f{Xs), 


Js{x) = 


j d'^X2d'^Xif{x 


- Xi)f{x2 


- x)tlj^{x2)[C"'^,T2f3^]tlj{x,), 


Ja{x) = 


j d'^Xid'^x^f^x 


- X3)fiXi 


- x)V'(Xi) [yf,CTiT2(3^]lp'^{Xs), 


Ja{x) = 


j d^X2d'^X4^f{x 


- X4)f{x2 


- x)^^J^{x2)[C''Yr2r^f3^]^P{x4 



(17) 

The object C = 27275 is the charge conjugation matrix. The matrices (3"^ = a/3/2A^(A = 
2,5,7) project onto the colour 3 channel with normalisation ti(l3"^l3^') = 35^^'. The cou- 
plings Qs and Qa specify the strength in the scalar and axial-vector diquark channels, respec- 
tively. We fix them by the empirical nucleon mass. The form factor in the diquark sector is 
assumed to be the same as in the meson sector. 

n the ladder approximation the quark-quark T-matrix near the pole can be parametrised 



as 



23 



24| 



T(pi,p2; h, k2) = Tip, k; q) ^ x\p. q)D\q)t{k, q) + q)D^'- {q)^ q) 



;i8) 



where q = Pi + P2 = ki + k2 is the total momentum of the quark-quark pair and p, k 
are the relative momenta within the diquarks. The momentum-dependent scalar x^i the 
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axial x'^ diquark vertex functions, their adjoint functions and x'^ and the corresponding 
)ropagators D^{q) and D'^'^{q) can be obtained from solving the Bethe-Salpeter equation 
3] for qq scattering matrix shown in Fig. 3. 

The scalar and the axial diquark masses Mds and M^a are obtained as the position of 
the pole of the corresponding T-matrix in the scalar and axial-vector diquark channels, 
respectively. 

With the diquark, the baryon can be constructed as a bound state of diquark-quark. We 
ignore the three-body irreducible graphs. Due to the separability of the two-body interaction 
in momentum-space, the relativistic Faddeev equation can be recast into an effective two- 
body Bethe-Salpeter type equation (pictorially shown also in Fig. 3) 



where the G~^{p,k\K) denotes the inverse of the full quark-diquark 4-point function and 
contains the sum of the disconnected part and the one-quark exchange interaction kernel. 
After projecting the effective diquark-quark Bethe-Salpeter equation in the colour singlet 
and isospin 1/2 channel, one finds 



where x^ and x^'^ stand for the Dirac structures of the diquark vertices introduced in Eq. (fTH|) . 
The dressed quark propagator S~^{p) is defined in Eq. (jH)). The numerical coefficient in kernel 
of the Faddeev equation (j2(Jj) comes from projecting the kernel to the physical baryon states 
(colour singlet and isospin half). We define the spectator quark momentum pq = r]P + p and 
the diquark momentum pd = {l—i])P—p, where P is the total momenta in the diquark-quark 
pair and rj G [0, 1] is the Mandelstam parameter which parametrises the relative momenta 
within the quark-diquark system^. The relative momentum of quarks in the diquarks vertices 
are defined aspi = p+k/2 — {l — 3ri)P/2 andp2 = —fc—p/2+(l—3r7)P/2, and the momentum 

^ Notice that observables do not depend on the Mandelstam parameters rj . 




(19) 



G 



1 




(20) 



12 



of the exchanged quark is g = —p — A; + (1 — 2ri)P. We refer to Appendix for details of 
numerical methods involved in solving the Faddeev equation ((201) • 

The parameters of the models Qni^c and A are fixed by meson properties i.e. the pion 
mass and decay in vacuum as described in section II. We select parameter set C given in 
Table 1. Now the only unknown parameters are diquark couplings Qs and ga and vector 
meson coupling g^. The scalar and axial- vector diquark couplings are taken as independent 
parameters. One can find a line in parameter space of diquark couplings in which a rea- 
sonable description of nucleon does exist, (see Fig. 11 in Ref. [IS]). In other words, the 
interaction is shared between the scalar and the axial-vector diquark and for small scalar 
diquark coupling gs one needs a dominant axial-vector diquark ga > gs and inverse. We 
fix these diquark couplings in vacuum to obtain a nucleon mass of 940 MeV. We choose 
gs = O.ligj, and ga = 0.23(7^ (the value of g,, for set C is given in Table 1) and corresponds 
to a scalar diquark with a mass M^s = 705 MeV and an axial diquark with a mass M^a = 660 
MeV in vacuum hm. Note also that it has been shown that the axial- vector channel is much 

n 

more important for the confining than non-confining parameter sets of our model [1^. All 
the results in the next section are given by the above-mentioned parameter set. The vector 
coupling g^ is the only free parameter which is left for the finite density calculation. 



IV. NUCLEAR MATTER 



In general, using a functional integration technique (similar to bosonization scheme pre- 
sented in section II) one can recast the Lagrangian in terms of hadron degrees of freedom 
and some auxiliary fields such as diquarks. The quark and the auxiliary fields can be inte- 
grated out and a chemical potential for the nucleons should be introduced. Having done that 
one can then directly apply conventional many-body techniques to the effective Lagrangian. 
Following Bentz and Thomas, we construct baryons out of the diquark-quark loop by solv- 
ing the Faddeev equations and solve baryonic matter in a mean-field approach^. Since we 
already have ignored meson loops, we do not keep diquark loops as welF. The vacuum 

^ We use a simple approximation that the nuclear matter expectation value of any operator, consists of its 
expectation value in the vacuum and an average over the nucleon Fermi-sea of correlated valence nucleon 

a. 

Diquarks are alike mesons boson and furthermore, our consideration is limited to a low-density region 
where diquarks do not condensate. 
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contribution of quark fields is kept at one- loop level. Therefore, the effective potential fl{p) 
regularised by subtracting its corresponding value at zero density can be written as 

Q{p) = fivac + ^^Med, (21) 

the contribution due to vacuum part i^vac is 

^vac = / 7^ In ,, ,JJ, + - -j^- (22) 



(27r)4 k^-M^ik) 2g^ Ag^ 
The medium part contains the nucleons and is defined 

/d^k 
j^^EN{k)nN{k), (23) 

where 7Ar = 4 is the spin-isospin degeneracy factor for the nucleon. The eriergy spectrum of 
a single nucleon has a simple form E^ik) = a/A;^ + + 3uf (for = 0) |5j] where Mj\f{'a) is 
the effective nucleon mass in medium obtained by solving the bound state Faddeev equation 
()19|) . The nucleon mass now is a complicated function of scalar field a. We define the baryon 
density p by 

/d^k 
(^^^(^)' (24) 

where the Fermi distribution function niv(fc) is at zero temperature the step function nN{k) = 
e{kF-k). 

Imposing the self-consistence condition ^ = leads to the following finite-density gap 
equation 

_ . f d^k M{k)p{k) f d^k MnOMm 

a = 2(7,7,2 / 77^-77 ' 9nlN / ^ ^ nN{k). (25) 



{2nyk^-M^{k) "^"'J {2nfEN{k) da 
This is a highly non-linear equation which is now coupled with the solutions of the Faddeev 
equations. Therefore, in order to solve the gap equation at a given density one needs to 
solve the Faddeev equation at the same time. In this way, the non-perturbative feature of 
the nucleon as a composite object is taken into account. This equation resembles the gap 
equation in the vacuum Eq. (fTUI) . In order to obtain the physical value of aJo, we require 
that = 0, for nuclear matter at rest this yields 

/d^k 
^^-pn(A;) = Qg^p. (26) 

Notice that the general form of the Faddeev equation ()19|) for the nucleon and the Bethe- 
Salpeter equation for the diquark in nuclear medium remains the same as in the vacuum. 
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FIG. 4: On the left the zero-momentum constituent quark mass M(0) and the (real, imaginary) 
part of the first pole of the dressed quark propagator (M/j, Mj) are displayed as a function of 
density p in the unit of the nuclear matter density po . We also show the effective nucleon Mj\f and 
the diquark masses Mfig,Mfia (for = 0) at finite density. On the right we show the fictitious 
scalar, axial diquark-quark threshold Mds + Mn and M^v+Mfi, respectively as a function of density 
p/po (the nucleon mass Mn is also shown). 

The only new input is the in-medium modified scalar mean field Eq. ()25|) which changes 
quark and diquark propagators appearing in Eq. (pUj) . At low density we do not have quark 
matter background, the Pauli-blocking is taken into account in constructing the composite 



nucleon and in the nuclear matter , lU | . 

In Fig. 4 (left panel) we show the effective quark mass obtained as the solution of the 
finite-density gap equation ()25p. the scalar diquark mass M^s and the axial-diquark mass 
Mda at finite density. The density dependent nucleon mass obtained from Faddeev Eq. 
is also shown. Notice that the nucleon mass does not follow from the addition of the diquark 
and constituent quark masses because of the relativistic Faddeev equation. 

The effective constituent quark mass in the nuclear matter decreases with a different 
slope compared to that in the quark matter background . This indicates that the many- 
body effect felt by a constituent quark in the nucleon is different from the quark matter. In 
Fig. 4 (left panel) we also show the real M^(p) and imaginary Mj{p) part of the first pole 
of the dressed quark propagator as a function of nuclear matter density. The deconfinement 
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occurs where Mj{p) = 0, this corresponds to p/po = 0.28. The real part of the first pole of 
the dressed quark propagator Mr increases with density up to the deconfinement point and 
then it decreases and approaches the zero-momentum constituent quark mass M(0) near 
the nuclear matter density po- In our model, we do not have a well-defined quark-diquark 
threshold. However, one may define a fictitious quark-diquark threshold as Mds + Mr (or 
Mda + Mfi), where Mr is the real part of the first pole of the dressed quark propagator. It 
is noted from Fig .4 (right panel) that the nucleon mass remains below the fictitious scalar 
diquark-quark threshold for all densities. However, it moves slightly above the fictitious 
axial diquark-quark threshold near the nuclear matter density po. Note that we show in 
Fig. 4 only the real part of nucleon mass. Above the diquark-quark threshold nucleon pole 
position is moved in the complex plane where the real part is taken as physical mass. Our 
numerical procedure (see Appendix) is valid when nucleon solution is not far above the 
fictitious threshold. It is natural to assume that at some point as we increase the density, 
nucleons decay to quarks and diquarks due to deconfinement. At the moment, we do not 
know yet the order of a possible deconfinement transition and its critical density. Based on 
our model approach, one can not address such questions quantitatively. 

The scalar and the axial-diquark masses in medium decrease in a very similar way. There- 
fore, their relative importance in the nucleon description remains almost the same at different 
densities. The nucleon mass compared with the diquarks and the quarks has a bigger de- 
creasing slope at very low density. However, it tends to saturate near the nuclear matter 
density. This effect is also observed for the diquarks, however, less pronounced. 

The stabilisation of the nucleon mass at about nuclear matter density and the dy- 
namical partial restoration of the chiral symmetry may closely be interconnected. This can 
be realized, since after chiral restoration of the constituent quark mass, the nucleon cannot 
feel the presence of the nuclear medium any more. The medium is incorporated via the 
quark propagator Eq. (jH)) through the finite-density gap equation p5|l into the diquark and 
consequently in the diquark-quark Bethe-Salpeter equations. The Pauli-blocking effect in 
the medium is not present for the quarks since in the medium we do not have a quark matter 
background. We employ the quarks only to construct the nucleons and assume that they are 
not resolved by the medium. The stabilisation of the nucleon mass in nuclear matter is the 
same as before. In the local NJL model where quark confinement is simulated by introduc- 
ing an infrared cutoff the stability of the nucleon mass is reached at higher density j^. In 
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FIG. 5: We show the density function of the scalar and axial diquark components within the 
nucleon at a given nuclear matter density p with respect to relative momentum between diquark 
and quark. All given components have spin 5=1/2 but different orbital angular momenta L as 
indicated in the plot. 

the context of the derivative scalar coupling model this phenomenon is related to dynamical 
screening of the effective coupling ji37|]. In our model, it is derived from the non-locality of 
the underlying theory. 

In order to find a better understanding of the internal structure of the nucleon in medium, 
we co„.t.uct the speetato. ,„o„,e„tun, d,stntatio„ fu„ct,o„ with. „.e.eon Q a. a 

function of nuclear matter density in terms of various diquark components in the nucleon 
rest frame. 



p{\p\,M^) = j dp,^^\\p\,p,,M^){D\p,))-'ij\\p\,p,,M^ 



N 



(27) 



where ip^^ and ip'^" are the Faddeev vertex functions obtained from Eq. (fT^ and p stands 
for the space component of relative momentum p. D^lpa) and D^^{pd) denote the diquark 
propagators. Although this definition is not unique, it provides some useful information 
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about the nature of compositeness of the nucleon and their response to the medium. Note 
that the nucleon vertex function contains information about the compositeness of the diquark 
through the diquark vertex function x^'^, see Eqs. fll8p2()j) . 

In the rest frame the mass of the nucleon and its total angular momentum are good 
quantum numbers. In the rest frame it is also possible to decompose the nucleon vertex 
function in terms of tri-spinor in different diquark channels each possessing definite orbital 
ang ular momentum and spin which allows a direct interpretation of different components 



|18l |22| . The results are shown in Fig. 5. It is observed that at finite density the s- 
wave in the scalar diquark channel is the dominant contribution to the nucleon ground 
state. As the nuclear matter density increases the scalar and the axial diquark masses 
decrease, consequently the diquark interaction couplings Qs and Qa grow. This effect can be 
realized from Fig. 5, since the maximum strength of the density function increases as we 
increase the baryonic density. It is interesting to notice that from Fig. 5 one can observe 
that the relative importance of the scalar and axial diquarks remains intact in medium (at 
least at low density). In order to estimate the change of the nucleon size in response to 
nuclear matter environment, we compute Prms = ((p^) ~ (p)^)^^^ at finite density. This 
leads to pRMs(MeV) = 232,182,146 for p/po = 0,0.22 and 0.40, respectively. Therefore, 
as baryonic density increases the nucleon size grows. An increase about 35% is found for 
density p/po ^ 0.40. 

Fig. 6 shows the binding energy per nucleon as a function of density. The equation 
of state when the internal structure of nucleon are taken into account in a presence of a 
vector field indeed saturates. We adjust the vector mean field coupling in order to have 
right binding energy per particle Eb/A = —16 MeV. This corresponds to p/po = 0.60 
when = 5f^/2.4. The fact that the equation of state curve cannot pass exactly through 
the empirical saturation point implies that our hybrid approach based on a mean field 
approximation in the medium is still a very crude approximation^. 

It is well known that in the chiral models the stabilitv of the nuclear matter depends 



very much on the dynamical chiral restoration 



45l |. The necessary condition for 



saturation of nuclear matter in any relativistic mean field theory is that nucleon attraction 



Notice that here we have only one free parameter in medium and it is not straightforward to fix two 
values via it. A similar study by Bentz and Thomas in the context of the local NJL model has shown 
that the binding energy saturates at too high density i5]. 
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FIG. 6: Energy per baryon as a function of p/ {po is nuclear matter density), for parameter set 
C, when = 5^/2.4. 



mediated from the a-meson exchange should decrease at high density. Therefore, if cx-meson 
mass decreases too rapidly with de nsity due to chiral restoration, it may work against the 
stabilization of the system 

mm 

However, as we have numerically proven, it is in 
principle possible to obtain a saturating nuclear matter equation of state, if the scalar field 
couples with the quarks instead of directly with the nucleons (see also [^). To this end, the 
nucleon is taken as a diquark-quark state which moves in self consistent scalar and vector 
fields coupling to the quarks. 



V. CONCLUSION 



In this paper we have investigated the nuclear matter in a non-local NJL model. We have 
worked out the equation of state of the nuclear matter in a hybrid approximation. First, we 
constructed the nucleon as a bound state of a diquark and a quark in the relativistic Faddeev 
approach and then built up nuclear matter by means of these composite nucleons. In this 
way, we incorporated the internal quark structure of the nucleon in a very simple framework. 
We showed that in this approximation the binding energy saturates. This is in contrast to 
the local NJL model, where in a similar framework, a stable normal nuclear matter does 
not exist unless one introduces a new extra parameter into the model either by including an 



19 

8-Fermi interaction term P| to induce the original coupling of the model density dependent 
or an infrared cutoff to simulate the confinement effect ^,|4^. Therefore, the long standing 
problem of matter stability in the NJL model can be resolved by introducing non-locality 
without invoking any adhoc new parameters. Note that the form factor Eq. (jlHj) does not 
introduce any new parameter except the cutoff A which is also present in the local NJL 
model. 

We have also studied the nucleon properties such as the modification of the nucleon 
mass and size in a nuclear matter medium. The nucleon mass in the medium decreases 
very rapidly but saturates near the nuclear matter density po- We obtained the nucleon 
wave function from the Faddeev equation and showed that the nucleon size significantly 
increases in the medium. This imphes a swollen nucleon in nuclear matter which has many 



391, |4 



41[. Our estimation of in-medium nucleon size should 



attendant consequences |3 
be taken more qualitatively since the quark density function given in Eq. ()27|) has not a 
unique definition. We found that the swelling of the nucleon is about 35% at about half 
nucleon matter density. One should note that in our model at nuclear matter saturation 
density p = O.Gpo the "confinement" mechanism is no longer at work, see Fig. 4 (left panel) 
and the nucleon becomes close to the fictitious diquark-quark threshold, see Fig. 4 (right 
panel). Moreover, since the nucleon mass saturates at high density, the nucleon will not 
expand forever^. Other studies have predicted different values^". For an example, it has 
been shown within the Friedberg-Lee nontopological-soliton bag model that the swelling of 
the nucleon is about 30% at normal nuclear matter density [ssl- In the Skyrmion picture 
nucleon swelling about 40% has been reported [4^. It has been also shown that in the 
quark-meson coupling model nucleon swelling about 25 — 35% in saturated nuclear matter 
can explain the observed depletion of the structure function in the medium Bjorken x region 



4l|. 



We have also investigated the role of diquarks within the nucleon in the nuclear matter 
medium. The scalar and axial diquark masses decrease with a very similar slope in the 
medium. Despite the fact that the nucleon mass decreases (and its size increases) signifi- 
cantly, the role of the scalar and the axial diquark in the nucleon description remains almost 



^ We could not obtain the Faddeev wave function at the higher density due to computation difficulties to 
verify this effect. 

^ One should be aware that nucleon size has not generally a unique definition among various models. 
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the same in the nuclear medium. 

In this paper, for simphcity we assumed that the auxially vector field to be local. The role 
of the vector meson in medium is still an open question and it deserves further investigation 
in our model as well. 
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APPENDIX: NUMERICAL METHOD 



There are some subtleties involved in evaluating loop integrals as of Eq. (fTT|) in Euclidean 
space. For simplicity this integral is evaluated at the timelike momentum q = (0,igo) in 
Euclidean space. For the confining parameter set, each quark propagator has a pair of 
complex conjugate poles. As we increase go; these poles in S{p±) (where p± = p ± |g) 
have a chance to cross the real axis and may produce an imaginary part in the meson (or 
diquark) propagator when a threshold for decay of a meson (or diquark) into qq- unphysical 
states opens^^, see Fig. 9. Therefore, the standard Wick rotation of the integration contour 
cannot be applied. We use the prescription proposed by Cutkosky et al. and successfully 
applied to such models in Refs. 13, 16, 4^. This amounts to a deformation of the integration 
contour (as indicated in Fig. 10) to ensure that the meson or diquark propagator does not 
develop an imaginary part, when the first pair of complex poles of quark propagator cross 
the real axis. At the pinch point, both the naive integral over Euclidean four-momentum and 
;he residue contribution diverge, although these divergences cancel to leave a finite result 



13 



421. 



The diquark-quark loops in the Faddeev equation (fT^ are computed in the rest frame of 
the nucleon P = {E, 0). First, we identify the singularities involved in the kernel of Faddeev 



The quark propagator has many set of quartets of complex poles and for the deconfinement case we have 
purely real poles in form of doublets. The external momentum going into the gg-loops (or gg-loops) does 
not change the structure of the quark propagator poles, however it changes their positions. In other words, 
if the relation Eq. (|14|l is satisfied with the parameters of the model, poles will always remain in form of 
quartets jl3| . 
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FIG. 7: We show the first complex conjugate poles of loop integrals in the p4 plane. On the left 
panel, the pole positions at qq = are shown. The open boxes and filled circles denote the poles of 
the p- and respectively. On the right panel, we show the deformed integration contour beyond 
the pinch point. 

equation (pUj) with respect to the integration momentum /cq- The singularities from quark 
propagator will occur at : 



pole 1 : 
pole 2 : 



ko = -r]E + ujI- ie, 
ko = -r]E - u;l + ^e, 



(A.l) 
(A.2) 



where tul = \/k^ + M^. The singularities from diquark propagator occurs: 



pole 3 : fco = (1 — r])E + cu'^ — ie, 

pole 4 : ko = {l-r])E -u'l + ie, 



(A.3) 
(A.4) 



where we define = \/k'^ + and is the the lower-energy pole in the diquark T- 
matrix (the mass of the lighter diquark) . The exchanged quark propagator is also singular 
at: 



pole 5 : ko = -po + 0;^+^ + (1 - 2?7)E - ie, 
pole 6 : ko = -po - cu^,^^ + (1 - 2r7)-E + ie. 



(A.5) 
(A.6) 



with uj^^^ = y/{k + py + M^. Note that all masses in the above equations are density 
dependent. We define 



a(p) = (1 - v)MRip) - vM.ip) r] e [0, 1]. 



(A.7) 
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It is easy to verify that in the complex k^, the poles 2 and 4 lie always left of fco = ct(p) and 
poles 1 and 3 lie always right of ko = a for E < Mji + M^- To avoid the singularities of 
the kernel, we perform a Wick rotation at a given density on the energy variables around 
ko = a{p), i. e. ko a{p) + iki (— oo < ki < +oo). Although a changes with density, 
it is in principle independent of the Mandelstam parameter rj. This provides an internal 
check on numerical procedure. In Refs. ^, Q|, the parameter rj is taken 1/2 from the 
outset. The location of the quark exchanged poles 5 and 6 depends on the po. One needs 
to bend the integration contour in such a way that the pole 5 and 6 to lie right and left of 
the path, respectively. To this end, we continue the momentum po to complex plane with 
Po = a{p) + ik2 then the poles 5 and 6 lie always right and left of the integration path 
ko = a{p) + iki. Above the diquark-quark threshold the solution of the effective Faddeev 
equation should be searched in complex P^-plane. In this case the complex pole coordinates 
p2 _ ^2 _ interpreted as a mass, E, and 7 is related to a corresponding width 

Our numerical method is valid when the location of the complex pole is close to the real 
axis and we concentrate only on computation of the mass. In practice, one may increase 
the stability of the numerical procedure by taking advantage of r] (in order to escape the 
above-mentioned poles far enough away). This implies that one can find aplateau within 



the range < 77 < 1 where the results do not depend on the parameter t] |]J 



22, M. The 



resulting complex coupled integral equations are solved by using the algorithm introduced 



by Oettel et al. 
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